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I.  Introduction 

The  turbulent  boundary  layer  at  constant  pressure  has  been  the  subject 
of  experimental  and  theoretical  investigations  for  many  years,  and  provides 
a well  documented  flow  for  the  assessment  of  experimental  techniques.  These 
investigations  have  shown  that  for  the  incompressible  turbulent  boundary  layer, 
the  turbulent  shearing  stress  can  be  measured  directly  or  can  be  calculated 
from  the  distribution  of  mean  velocity  with  the  aid  of  well-established  simi- 
larity laws.  For  compressible  flows,  however,  measurements  of  the  Reynolds 
stresses  are  rare,  and  interpretation  of  the  results  is  difficult.  Recent 
measurements  by  Johnson  and  Rose  (1973),  Yanta  and  Lee  (1974),  and  by  Abbiss 
(1976)  have  used  laser-Doppler  velocimetry  techniques  to  make  direct  measure- 
ments of  the  Reynolds  stress  in  turbulent  boundary  layers  with  free-stream 
Mach  numbers  in  the  range  1.3  to  3.0.  However,  a serious  anomaly  is  exhibited 
by  these  measurements  in  that  the  maximum  value  of  - pu  'v  ' occurs  much  further 
from  the  wall  than  is  reasonable  for  flow  at  constant  pressure.  This  anomaly 
has  been  discussed  by  Sandborn  (1974),  who  supports  the  conjecture  by  some  of 
the  authors  cited  that  density  fluctuations  may  contribute  substantially  to 
the  turbulent  stresses  near  the  wall.  This  conjecture  is  in  direct  opposition 
to  the  conclusion  by  Morkovin  (1961)  that  effects  of  density  fluctuations 
should  be  small  compared  to  effects  of  variations  in  mean  density  for  Mach 
numbers  up  to  4 or  5. 

The  purpose  of  the  present  experiments  is  to  obtain  redundant  data  over 

a substantial  range  of  Mach  numbers  (Mg  “ 0.1  to  2.2),  in  an  effort  to  resolve 

the  anomaly  in  turbulent  shearing  stress.  The  low-speed  experiments  were 

ie 

performed  in  the  boundary  layer  on  a flat  plate  model  in  the  GALC1T  Merrill 

* 

Graduate  Aeronautical  Laboratories  of  the  California  Institute  of  Technology, 
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wind  tunnel  at  a free-stream  Mach  number  0.1.  The  high-speed  experi- 
ments were  performed  in  the  ceiling  boundary  layer  of  the  J.P,L.+  20-inch 
wind  tunnel  at  Mach  numbers  0.6  s 2.2. 

The  flow  was  documented  by  conventional  means  using  a Pitot  tube, 
which  was  traversed  through  the  boundary  layer,  to  measure  the  mean  flow. 

In  addition,  surface-friction  measurements  were  made  using  both  a floating- 
element  balance  and  a Preston  tube.  The  mean-flow  scaling  suggested  by 
Van  Driest  (1955)  was  applied  to  the  data,  to  test  the  adequacy  of  a single 
similarity  formulation  for  both  compressible  and  incompressible  flow.  The 
shearing  stress  was  computed  from  the  mean  flow  as  part  of  the  analysis. 

The  details  of  these  experiments  have  been  discussed  in  Part  I of  this 
report  by  Collins,  Coles  and  Hicks  (1978). 

In  the  present  report,  the  instrumentation,  data  acquisition  and 
analysis  for  laser-Doppler  velocity  measurements  in  the  same  flows  are 

1 2 ,2  — j y 

described.  Measurements  of  u,  v,  u , v and  the  Reynolds  stress  - pu  v 
were  made  over  the  full  Ifech, number  range.  These  measurements  are  discussed 
and  compared  to  the  results  for  the  mean  flow  presented  in  Part  I of  this 
work . 

II.  Measurement  Considerations  at  High  Speeds 

The  use  of  laser-Doppler  velocimetry  techniques  in  high-speed  flows 
unavoidably  results  in  very  high  frequency  signals.  Figure  1 is  a plot  of 
the  Doppler  frequency 


= 2 ain®/2  u , 
u \ J- 


Jet  Propulsion  Laboratory 


(1) 
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as  a function  of  the  flow  velocity  u±  with  the  beam  separation  angle  6 
as  a parameter.  For  these  computations!  the  wavelength  has  been  assumed 
to  be  \ ■ 0.5145  fitn,  corresponding  to  an  argon-ion  laser. 

A reduction  in  the  beam  separation  angle  6,  in  order  to  achieve  a 
reduction  in  the  signal  frequency!  results  in  a reduction  in  the  number 
of  fringes  within  the  measurement  ellipsoid,  and  a consequent  reduction 
in  the  measurement  accuracy.  The  interesting  beam  pair  produces  a set  of 
virtual  fringes  whose  spacing  is  given  by 


1 

2 sin3/2 


(2) 


Because  the  (TEM^  Gaussian  beam  radius  at  the  focus  is  given  by 


v 

o 


15 

TTW 


I 


(3) 


where  R is  the  beam  wavefront  radius  of  curvature  after  the  focusing  lens 

2 

(R  r>J  f ■ lens  focal  length),  and  u is  the  L/e  intensity  envelope  radius 
at  the  focusing  lens,  the  resulting  number  of  fringes  is  then  given  by 
(e.g.,  Dimotakia  (1976)) 

w 

- 4 tan6/2  . (4) 

The  number  of  fringes,  1b  plotted  in  Figure  2 as  a function  of  the 
focal  volume  diameter  2 wq  with  the  angle  6 as  a parameter  for  a wavelength 
\ = 0.5145  urn.  As  can  be  seen  from  Figure  2,  spatial  resolution  consid- 
erations place  a lower  bound  on  the  angle  6,  presenting  a criterion  which 
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conflicts  with  the  requirement  to  keep  the  resulting  Doppler  frequencies 
small. 

Other  difficulties  arise  if  the  Doppler  frequencies  are  permitted 
to  be  large.  In  addition  to  the  problems  of  handling  high  frequencies 
per  se,  the  relative  accuracy  5u/u  with  which  the  velocity  of  a single 
particle  can  be  measured  is  never  less  than  the  ratio  of  the  processor 
clock  period,  t , to  the  total  flight  time  At  that  is  used  for  the  mea- 
surement, i.e.. 


6u 

u 


(5) 


By  way  of  example,  a processor  with  a clock  frequency  of  100  MHz 
(t^  = 10  sec)  can  measure  the  velocity  of  a Mach  ~ 2 particle  in  air 
(u^  ~ 530  m/sec)  in  a 800  pm  diameter  focal  volume  to  about  ~ 0.7%. 

In  addition  to  these  considerations,  which  would  limit  the  accuracy 
even  if  the  measurement  environment  was  noise  free,  other  factors  become 
important  at  high  velocities  by  decreasing  the  signal- to-noise  ratio  of 
a single  reading.  As  the  flow  velocities  become  higher,  smaller  particles 
are  required  to  minimize  problems  of  particle  lag.  This  results  in  rapidly 
decreasing  scattering  cross-sections  (Born  and  Wolf  (1975))  and  consequent 
signal  intensities.  In  addition,  as  the  velocity  increases,  the  scattering 
particle  spends  less  time  in  the  focal  volume.  As  a result,  even  though 
the  number  of  scattered  photons  per  second  remains  constant,  proportional 
to  the  scattering  cross-section,  fewer  total  photons  are  collected  by  the 
detecting  optics.  An  equivalent  way  of  analyzing  this  problem  is  to  consider 
the  noise  bandwidth  that  must  be  utilized  as  the  velocities,  and  resulting 
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Doppler  frequencies.  Increase,  For  a shot-noise  limited  detection  process, 
the  noise  power  increases  proportionally  to  the  bandwidth  that  is  accepted 
by  the  processing  electronics.  Consequently,  using  the  same  particles, 
the  s ignal-to-noise  ratio  will  decrease  as  the  velocity  increases. 

III.  Scattering  Particles 

A.  Particle  Lag 

The  extent  to  which  the  measurement  of  particle  velocity  represents 
the  fluid  velocity  is,  of  course,  a separate  issue.  If  the  Lagrangian 
fluctuation  frequencies  in  the  fluid  are  small  compared  to  the  reciprocal 
of  a characteristic  particle  response  time,  t , the  particles  are  expected 
to  follow  those  fluctuations. 

For  small  solid  particles  in  gas  flows,  the  equation  of  motion  in 
a Lagrangian  frame  can  be  approximated  by 


d 

dt  Up 


(6) 


If  the  particle  Reynolds  number 


Re 

P 


(7) 


is  small,  and  the  flow  can  be  described  by  the  continuum  equations,  we 
have  the  classical  Stokes  flow  for  a spherical  particle  for  which  the 
drag  coefficient  is  = 24/Re^.  The  particle  response  time  is  then  given 
by 


(8) 
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which  is  independent  of  the  fluid  density.  A plot  of  the  particle  response 
time  as  a function  of  the  particle  diameter  is  given  in  Figure  3 for  di- 
butyl phthalate  aerosol  particles  in  air. 

In  order  to  estimate  the  fluctuation  frequencies  to  which  a particle 
will  be  subjected,  it  is  useful  to  examine  the  Eulerian  fluid  motions  in 
the  turbulent  boundary  layer.  From  the  measurements  of  Klebanoff  (1954) 
in  a boundary  layer  at  constant  pressure,  and  from  the  measurements  of  Perry 
and  Abell  (1975 j 1977)  in  pipes,  the  behavior  of  the  Eulerian  velocity 
fluctuation  spectrum  can  be  inferred.  These  results  give  spectra  which  can 
be  described  by 


u'2(y) 


tt6  f r u)5  1 

u(y)  °U(y)J 


1 


for  the  wake  region,  and 


S (tu.y) 

uu 

U (y) 


(9) 


(10) 


in  the  inner  portion  of  the  boundary  layer,  but  outside  the  sublayer. 
These  spectra  then  result  in  expected  frequencies  given  by 


(u)>6  * const. 


(ID 


in  the  wake  region,  and 


(cu)  = const,  u(y) 

y y 


(12) 
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in  the  boundary  layer  interior.  Consequently,  particles  above  the  sublayer 
should  track  the  flow  well,  provided  that 


si  , 

y 


(13) 


or 


Lls  < JL  r y/6  1 . 

5 Me  U(y)/uJ 


(14) 


Substituting  Equation  (8)  for  7^,  we  obtain  an  upper  bound  on  the  particle 
diameter  for  the  particles  to  track  the  flow. 


d 

P 


< 


d 

Pinax 


«s 


(15) 


for  di-butyl  phthalate  aerosol  in  air  at  2S°C. 

A prediction  of  Equation  (12)  is  that  the  turbulent  wavenumbers  k rv  w/u 
increase  as  the  wall  is  approached,  . 


k 


const . 

y 


(16) 


This  result  is  confirmed  qualitatively  by  the  photograph  in  Figure  4 in 
which  the  turbulent  boundary  layer  on  the  surface  of  a water  channel  has 
been  made  visible  using  a thick  suspension  of  aluminum  flakes  which  pref- 
erentially align  with  the  local  strain  field.  The  photograph  was  taken  in 
the  GALCIT  low-speed  water  channel  by  Brian  Cantwell  and  is  reproduced  here 
with  his  permission.  The  Reynolds  number  for  this  photograph,  based  on  x. 
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was  106.  The  monotonic  decrease  in  scale  size  as  the  wall  is  approached 
is  quite  evident.  These  observations,  and  the  results  of  Equation  (12)  are 
consistent  with  the  usual  mixing-length  hypotheses  made  for  the  outer  layer 
by  Prandtl  (1942),  and  are  consistent  with  the  hypothesis  for  the  region 
closer  to  the  wall,  but  outside  the  viscous  layer  by  von  Karman  (1930). 

In  the  interior  of  the  boundary  layer,  the  quantity  in  the  brackets 
In  Equation  (14)  is  bounded  by  the  relation 


(17) 


where  u is  the  friction  velocity  defined  as 

T 


(18) 


The  ratio  u has  been  shown  by  Coles  (1968)  to  be  a function  of 

e 

Re0  for  incompressible  flows  with  zero  pressure  gradient.  The  square  root 
of  this  quantity  is  plotted  in  Figure  5,  for  the  range  of  Reg  covered  by 
the  present  measurements  (8  x 10^  < < 4 X 10  ) . The  data  for  this  figure 

are  taken  from  Table  4 of  Collins,  Coles  and  Hicks  (1978) . These  data  are 
compared  in  Figure  5 to  a straight  line  fit  given  by  the  emperical  relation 


f(*ee) 


5.01  Reg 


-0.435 


(19) 


In  writing  the  lower  bound  in  Equation  (17),  effects  of  compressibility 
have  been  ignored.  Using  the  numerical  result  of  Equations  (15)  and  (19),  we 
then  have 
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d 

P 


< d 


max 


^(f)% 

Mg*  (pm) 
24.2 


far  from  the  sublayer 


near  but  outside  the  sublayer  . 


(20) 


By  way  of  example,  a flow  of  ~ 1 requires  particles  less  than  1.5  pm 
in  diameter  to  track  the  fluctuations  at  y/6  ~ 0.1. 

B.  Techniques  for  Seeding  the  Flow 

The  particle  generator  for  the  present  experiments  used  a Laskin 
nozzle  type  construction  to  generate  ,a  poly-dispersed  aerosol  of  di-butyl 
phthalate.  Particles  greater  than  1 pm  in  diameter  were  effectively 
removed  by  an  impact  plate  Incorporated  in  the  design.  The  resulting 
particle  size  distribution,  measured  using  both  a cascade  impactor  and 
a multi-channel  particle  analyzer,  is  shown  in  Figure  6.  These  results 
indicate  little  contribution  to  the  number  density  for  sizes  greater  than 
1 pm.  While  little  is  known  about  the  distribution  below  0.4  pm,  it  was 
possible  to  ascertain  that  no  particles  with  a diameter  less  than  0.1  pm 
were  present. 

IWo  different  methods  were  used  to  introduce  particles  into  the  flow 
depending  on  the  flow  facility  that  was  used.  For  the  high-speed  measurements 
(0.6  ~ Mg~  2,2),  which  were  conducted  in  the  JPL  20- inch  wind  tunnel,  it  was 
necessary  to  seed  the  flow  by  introducing  the  particles  into  the  settling 
chamber  between  the  last  turbulence  screen  and  the  contraction  section. 

The  particles  were  Introduced  through  a tube  which  protruded  5 cm  vertically 
into  the  settling  chamber  from  above.  Holes  were  drilled  into  the  tube  along 
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the  stagnation  line  as  well  as  the  rear  side  at  ± 30°  with  respect  to  the 
flow  vector.  Using  this  technique,  data  rates  in  excess  of  50,000  samples 
per  second  were  achieved.  Preston-tube  measurements  obtained  with  and 
without  particles  indicate  no  detectable  influence  of  the  introduction  of 

r 

particles  on  the  surface  friction  and  therefore,  probably,  on  the  mean  flow. 

For  the  low-speed  measurements  in  the  Merrill  wind  tunnel  (1^=“0.1), 
although  some  naturally  occurring  particles  were  present,  the  flow  was  seeded 
using  the  same  aerosol  generator  in  order  to  increase  the  data  rate  and  to 
control  the  scattering  particle  size  distribution.  For  these  measurements 
the  particles  were  Injected  through  a tube  spanning  the  test  section  located 
downstream  of  the  boundary- layer  plate,  providing  a uniform  distribution  of 
particles  throughout  the  test  section.  In  practice,  the  adjustment  of  the 
rate  at  which  particles  were  introduced  proved  difficult  because  the  charac- 
teristic time  was  of  the  order  of  several  minutes  for  equilibrium  with 
particle  losses.  As  a consequence,  the  particle  number  density  changed 
during  the  time  required  ,for  a traverse  through  the  boundary  layer.  Using 
this  seeding  technique,  data  rates  in  excess  of  2,000  samples/sec  were 
achieved  at  low  speeds . 

C.  The  Distribution  of  Particles 

The  particle  number  density  was  measured  by  counting  the  average  fre- 
quency of  occurrence  of  valid  data  from  the  processing  electronics.  These 
data  are  a measure  of  the  mean  particle  flux  (nu).  Assuming  that  correlations 
between  the  scattering  particle  number  density  and  the  streanwise  velocity 
are  small,  then 


(21) 
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gives  the  mean  value  of  the  particle  number  density. 

It  is  an  unfortunate  consequence  of  this  measurement  technique  that 
the  measured  value  of  the  particle  flux  (nu)  decreases  with  the  signal-to- 
noiae  ratio  of  the  experiment.  The  most  common  cause  of  a decrease  in  the 
signal-to-noise  ratio  in  the  high-speed  flow  experiments  has  been  attributed 
to  the  random  occurrence  of  oil  filaments  on  the  windows  through  which  the 
observations  were  made.  Imperfections  in  the  lucite  windows  in  the  Merrill 
wind  tunnel  had  a similar  effect  on  the  measured  particle  flux.  In  addition, 
a more  serious  problem  arose  in  the  low-speed  measurements,  where  the  particle 
flux  was  a slowly-varying  function  of  time. 

The  resulting  distribution  of  particle  number  density  through  the 
boundary  layer  is  shown  in  Figures  7,  8 and  9,  in  which  the  particle  number 
density  has  been  normalized  by  the  local  fluid  density  in  order  to  account 
for  the  effects  of  compressibility.  These  data  are  further  normalized  by 
their  maximum  value  in  the  boundary  layer  to  eliminate  the  effects  of  differing 
injection  rates. 

The  loss  mechanisms  of  the  scattering  particles  from  the  boundary  layer 
are  important  in  understanding  the  resulting  laser-Doppler  observations.  For 
the  high-speed  experiments,  the  particle  distribution  was  initially  both 
spatially  and  temporally  approximately  uniform  over  a restricted  volume  of 
fluid  in  the  settling  chamber.  The  particle  distribution  in  the  outer  portion 
of  the  boundary  layer  is  then  determined  by  the  diffusion  of  the  passive 
scalar  into  the  outer  flew.  In  addition  to  diffusion,  the  entrainment  of 
fluid  into  the  boundary  layer  from  the  free  stream  by  the  mechanism  of  large- 
scale  eddy  engulfment  contributes  substantially  to  the  structure  of  the  particle 
number  density  profiles.  Because  the  entrained  fluid  contains  no  particles. 
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the  measured  distributions  shown  in  Figures  7 and  B exhibit  a rapid  decay 
away  from  the  wall,  with  the  particle  number  density,  falling  to  zero  in 
the  outer  portion  of  the  boundary  layer.  This  entrained  fluid  is  under- 
represented in  the  Doppler  measurements  because  of  the  loss  of  statistical 
weight . 

As  the  wall  is  approached,  the  particle  number  density,  normalized 
by  the  local  fluid  density,  decays  more  rapidly  than  the  velocity  as  a result 
of  the  loss  of  particles  to  the  wall.  This  loss  to  the  wall  is  evident  from 
the  observation  of  an  oil  film  on  the  test  surface. 

IV . Optics  and  Measurement  Geometry 

A.  Position  Measurement 

The  present  experiments  utilized  the  laser-Doppler  velocimeter  in  the 
single  particle,  dual  forward-scatter  mode.  The  instrument  was  mounted  on 
a vibration  isolated,  two-axis  traverse  which  permitted  surveys  up  to  75  cm 
in  the  streamwise  direction,  and  20  cm  in  the  direction  normal  to  the  test 
plate.  For  measurements  in  the  JPL  20-inch  wind  tunnel,  the  instrument  was 
mounted  on  top  of  the  wind  tunnel  test  section,  as  depicted  in  Figure  10. 

For  measurements  in  the  GALCIT  Merrill  wind  tunnel,  the  traverse  was  suspended 
from  the  ceiling  of  the  laboratory,  and  was  not  fixed  relative  to  the  test 
plate . 

The  vertical  position  was  sensed  by  means  of  a helipot  mounted  on  the 
vertical  traverse  and  was  repeatable  to  better  than  0.005  cm.  The  origin  for 
the  vertical  position  measurements  was  estimated  by  placing  a razor  blade  of 
known  width  normal  to  the  test  surface  at  the  location  of  the  focal  volume. 

The  origin  was  then  measured  by  determining' the  half  intensity  of  transmission 
as  a function  of  position  of  one  of  the  beams. 
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In  practice,  establishing  the  origin,  to  within  the  accuracy  warranted 
by  the  data,  by  direct  measurement  proved  difficult  for  a variety  of  reasons. 
These  inherent  difficulties  in  measuring  the  position  of  the  wall,  relative 
to  the  LDV  traversing  mechanism  necessitated  establishing  the  origin  in  y 
using  the  data  for  the  mean-velocity  profile.  This  was  accomplished  by  a 
least-squares  fit  of  the  mean-velocity  data  to  the  law  of  the  wall  and  the 
law  of  the  wake  In  the  form 


+ 1 , + . , 211  , 2 v 
u = — In  y +c  + sin  Y , 
k H 

where 

u+  « u/ut  , 

+ UT 

y ■ (y  - yD)  “ * 

Y " *76  ■ 


(22) 


(23) 

(24) 

(25) 


and  where  y is  the  small  offset  of  the  position  of  the  wall  required  to 
o 

establish  the  proper  origin. 

For  this  purpose,  a least-squares  analysis  was  performed  to  determine 
the  parameters  , II  and  yQ  wi'th  k — 0.41  and  c “ 5.0  given,  and  6 was 
computed  as  a constraint  imposed  by  the  local  friction  law. 


+ 

u 

e 


—In  6+  + c 
K 


(26) 


The  technique  for  this  analysis  is  outlined  in  Part  I of  this  work  by 
Collins,  Coles  and  Hicks  (1978)  in  conjunction  with  the  mean-flow  data 
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presented  there.  The  results  of  this  procedure  lead  to  corrections  yQ 
of  the  order  of  the  focal  volume  diameter.  In  particular,  the  maximum 
correction  was  0.18  cm  for  one  case  of  the  high-speed  data  and  0.075  cm 
for  the  low-speed  flow  measurements. 

B.  Transmitting  Optics 

The  light  source  for  this  instrument  was  a Coherent  Radiation  Model 
52B,  4-watt  argon-ion  laser.  The  laser  was  operated  single  line  at 
0.5145  pm  and  used  an  oven  stabilized  etalon  to  provide  a single-mode  beam. 

It  was  mounted  on  the  tubular  support  structure  for  the  optics  as  shown  in 

l 

Figure  10.  This  configuration  ensures  that  the  focal  volume  remains  unchanged 
as  the  assembly  is  traversed  and  also  that  the  relative  alignment  between  the 
transmitting  and  receiving  optics  is  preserved  to  within  a few  percent  of  the 
focal-volume  diameter. 

The  laser  output  beam  was  directed  into  the  vertical  support  tube  for 
the  transmitting  optics  by  means  of  two  separate  mirrors  in  order  to  preserve 
the  proper  polarization  of  the  beam.  The  transmitting  optics  cell  contains 
a fixed  mirror  to  direct  the  beam  into  the  horizontal  plane  maintaining  a 
vertical  polarization  vector.  The  beam  is  then  split  into  two  beams  of 
unequal  intensity  by  reflection  from  the  front  and  rear  surfaces  of  the  first 
of  two  multiply-coated  beam  splitters.  The  resulting  beam  pair  is  then  trans- 
mitted through  a second  multiply-coated  beam  splitter  which  divides  the  stronger 
beam  into  two  beams  of  equal  intensity.  This  arrangement  results  in  three 
parallel  beams  of  approximately  equal  intensity  which  form  a right  isosceles 
triangle  whose  base  is  nominally  parallel  to  the  test  surface.  The  three 
output  beams  are  then  focused  in  the  center  of  the  tunnel  by  a 1.0  meter  focal 
length,  antireflection  coated  lens,  to  a common  focal  volume  ~ 0.8  mm  in 
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diameter.  The  optical  axis  intercepts  the  test  plate  from  below  at  an 
angle  of  ~ 1:100  in  order  to  permit  measurement  of  the  boundary  layer 
close  to  the  surface  of  the  polished  test  plate. 

The  resulting  focal  volume  contains  three  independent  sets  of  virtual 
fringe  planes  correspondingly  perpendicular  to  the  u,  (u_  + v) !*f  2 and 
(u  - v) /*J  2 velocity  vectors.  The  geometry  of  the  three  beam  pairs,  and 
actual  photographs  of  the  fringes  formed  by  each  beam  pair,  are  shown  in 
Figure  11.  By  selectively  blocking  one  of  the  three  beams,  any  one  of  the 
three  velocity  components  can  be  examined  without  a change  in  the  focal- 
volume  geometry.  This  method  we  consider  superior  to  the  more  common  one 
of  rotating  a single  beam  pair  in  the  desired  orientation  which,  in  our 
experience,  cannot  easily  be  done  with  sufficient  reproducibility  and 
without  slight  displacements  of  the  focal  volume. 

C.  Geometry  of  the  Intersection  Volume 

The  three  intersecting  beams,  taken  in  pairs,  define  three  planes, 
along  which  the  measurements  were  taken.  Even  though  considerable  care 
was  taken  to  orient  these  planes  at  0°  and  ± 45°  with  respect  to  the  wall, 
the  small  deviations  that  were  measured  were  taken  into  account  in  the  data 
analysis . 

The  subtended  angles  between  the  beam  pairs  were  measured  by  two 
different  methods  on  two  different  occasions  and  found  to  be  within  ~ 0.3% 
of  the  values: 

u measurement | =*  0.02112  rad 

y~-  (u  + v)  measurement 

- 1 ■ (ij  - v)  measurement 
*/2~ 


= 0.01490  rad  (27) 

= 0.01424  rad. 
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Thus,  ^ “ 8 define  a nearly  isosceles  triangle  as  illustrated 

in  Figure  11, 

As  a consequence  of  the  requirement  to  isolate  the  traversing  mechanism 
from  vibrations,  the  relative  orientation  between  the  hypotenuse  of  the 
triangle  and  the  wall  (angle  0,  Figure  12)  had  to  be  determined  by  a separate 
measurement  for  each  profile.  This  was  accomplished  by  permanently  mounting 
a second  razor  blade  to  one  side  of  the  test  plate  with  its  edge  aligned  with 
the  surface  of  the  plate.  The  relative  orientation  for  the  optics  could  then 
be  determined  by  measuring  the  offset  between  the  two  beams  forming  the  triangle 
hypotenuse.  This  angle  0 is  small  in  every  case  but  must  be  accurately  known 
to  estimate  any  quantities  that  involve  v,  the  velocity  component  normal  to 
the  wall. 

Using  the  geometry  as  defined  in  Figure  12,  we  have  from  the  law  of 


cosines , 


and 


cos  a = 


COS  P = 


cos  y = 


822  * - °12 
2 S2®3 

®32  + *1  - 622 

2 9^ 

2 2 2 
81  ^ -83 

2 V2 


(28) 


Assuming  now  that  the  velocity  vector  in  two  dimensions  can  be  defined 
by  its  components  in  the  three  directions  defined  in  Figure  12,  we  have 
that  the  three  measured  components  are  given  by 
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U1  = bllu  + b12V 


u2  = b21U  + b22V  > 


u3  " b31U  + b32V  y 


(29) 


where 


b^  = cos0 


b21  = C0S^Y  + ^ 


b31  “ cos<>^  " & 


b12  = sin0 


b22  = sin^Y  + & r 


b32  - -sin(3  - 0) 


(30) 


Equations  (29) form  a redundant  set,  and  car  be  used  to  check  the  validity 
of  the  data.  In  particular,  since  sin(P  + y)  = sinor,  we  must  have 


(u^>  sinor  = (u2>  sing  + <u^>  siny 


(31) 


Recall  from  Equation  (1)  that 


2 sin0./2  9. 

\ = Ui  “*  T Ui  ’ 


(32) 


so  that  the  consistency  Equation  (31)  becomes 


> - <vR  > + <v„  > 


(33) 


Alternatively,  Equations  (29)  can  be  solved  in  pairs  to  obtain  two  in- 
dependent estimates  of  (u)  and  (v)  at  each  location. 
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312\ 

b21 


and 
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31 


12 


32' 


/ 

X<vnV 


'<“l>' 


^u3'5' 


Inverting  Equations  (34)  and  (35),  we  then  have 

<ux)  sin(v  + 0)  - <u2>  sin0  "| 


S>  = 

<VI>  = 


siny 


<u2)  cos (0)  - (Uj)  cos (y  + 0) 
sinY 


and 


<UH> 


<VII>  = 


(ux)  sin(g  - ©)  + <u3)  sin0 
sing 

<u3)  coa (3  - 0)  - (u3>  COS0 
sing 


(34) 


(35) 


(36) 


(37) 


The  final  valueb  for  (u)  and  (v)  were  then  estimated  by  averaging  the  two 
independent  estimates, 


(u>  ~ J + (uxi^ 

<v>  - | ((V,)  + <»„>) 


(38) 


Using  Equations  (29),  we  can  also 
of  the  measured  velocity  components,  to 


re late 


the  mean-square  fluctuations 
(v  f^)  and  (u  fv  ')  . In  ^ 
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particular,  we  have 


<u'l2>  = bl;L2  <u/2>  + b122  <v'2>  +2  bxlb12 
<u'22)  - b212  <u'2>  + b222  (v'2>  + 2 b21b22 

<U32>=b312  <U'2>+b322  <v  ^ + 2 b31b32 


<u  'v  ') 


<u  v '>  y 


<u'v'> 


j 


(39) 


Note  that  while  Equations  (29)  are  redundant.  Equations  (39)  are  not, 
except  for  special  choices  of  g,  y and  0.  Inverting  Equations  (39), 
we  then  have  . 

<u'2>  = Cu  <u'i2>  + C12  <u/22>  + C13  <U,32> 

<v'2>  = C21  <u-)  + C22  (u'22)  + C23  (u^2)  " ‘ 

<»'»'>  "C31  <u  i ^ > + C32  <u'22>+C33  <u/32>^  * 


,(40) 


where  the  matrix  is  given  by. 


C11  = s^n(P  + Y)  sil*(g  - 0)  sin(y  + 0) /A 
C^2  = -sing  sin(g  - 0)  sin0/A 


C^3  = siny  sin(y  + 0)  sin0/A 


C21  = -sin(g  + y)  cos(g  - 0)  cos (y  + 0) / A 


C22  = sing  cos(g  - 0)  COS0/& 


C23  = siny  cos(y  + 0)  cos0/fi 


(41) 
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C31  " “2  sin^  + Y)  Sin(p  - Y * 2 Q)fb 
C32  = \ sinP  sin^  - 20)/i  ► 

C33  = slnY  sin^Y  + 2 0)/A  J ^ , 

and  where 

A = sin|3  • sinv  • sin(|3  + v)  . (42) 

D,  Receiving  Optics 

The  optical  axis  for  the  receiving  optics  is  aligned  at  an  angle  of 
~ 7.5°  with  respect  to  the  test-plate  surface,  thus  permitting  a full  view 
of  the  focal  volume  throughout  the  boundary  layer.  This  inclination  of 
both  the  transmitting  and  the  receiving  optics  with  respect  to  the  test 
surface  results  in  a reduction  of  the  spanwise  extent  of  the  focal  volume. 

The  collecting  lenses  consist  of  a pair  of  15.24  cm  diameter,  f/5  achromats 
which  have  been  anti-reflection  coated  and  are  mounted  back-to-back  in  a 
common  cell.  This  configuration  provides  a 1:1  imaging  system  that  was 
measured  to  be  approximately  twice  diffraction  limited.  The  collected 
light  cone  is  directed  by  a 10  cm  diameter  multi-layer  dielectric  coated 
mirror  into  a 0,05  cm  diameter  pinhole  which  spatially  filters  the  collected 
light  that  is  accepted  by  the  photomuLtiplier  assembly.  The  photomultiplier 
assembly  consists  of  a focusing  mount  with  a 20  mm  focal  length,  f/5  coated 
lens,  a 30A  wide  0.5145  pm  filter  and  magnetic  lensing  to  restrict  the  photo- 
cathode of  an  EMI  9813B  photomultiplier  tube.  The  photomultiplier  tube  output 
is  coupled  directly  to  a low-noise  fast  rise-time  preamplifier,  constructed 
inside  the  tube  housing  at  the  base  of  the  photomultiplier  tube.  The  photo- 
detection system  was  designed  with  a frequency  response  in  excess  of  50  MHz. 
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The  effective  focal -volume  diameter  is  defined  by  the  pinhole  in  the 
receiving  optics  and  is  therefore  equal  to  ~ 0.05  cm. 

Oscilloscope  traces  of  raw  signals  from  the  output  of  the  preamplifier 
at  the  base  of  the  photomultiplier  are  included  in  Figure  13.  These  photo- 
graphs are  multiple  exposures  taken  of  particle  bursts  at  the  highest 

4 

recorded  velocities  of  t^~2.2  (Re^  ~ 4 X 10  ) at  the  edge  of  the  boundary 
layer.  The  frequency  of  the  u-component  signal  is  about  26  MHz.  The  fre- 
quency of  the  (u  + v)-component  signal  is  almost  exactly  one-half  the  u- 
component  frequency.  This  feature,  which  will  be  discussed  later,  is  exploited 
by  the  processing  scheme.  Typical  signal- to- noise  ratios  of  15-20  db  were 
measured  at  these  velocities. 

V . Signal  Processing  and  Data  Acquisition 

The  output  of  the  photomultiplier,  coupled  through  the  integral  pre- 
amplifier, was  monitored  on  a HP  8553B/8552B  wave  analyzer  which  was  also 
used  to  align  the  optics  and  tune  the  entire  system  by  maximizing  the  signal- 
to- noise  ratio  of  the  Doppler  signal.  This  signal  was  filtered  by  a two-pole 
low-pass  filter  in  order  to  eliminate  the  substantial  shot  noise  above  the 
Doppler  frequencies.  The  pedestal,  resulting  from  the  dual-scatter  optics, 
is  not  removed  by  this  processing  scheme. 

A.  The  LDV  Processor 

The  filtered  output  signal  is  fed  into  the  LDV  processor  (Dimotakis  and 
Lang  (1977))  depicted  in  block  diagram  in  Figure  14.  The  various  components 
of  this  processor  are  described  in  the  following  paragraphs. 

1 . The  Analog  Input  Processor 

The  Analog  Input  Processor  defines  a reference  level  Vq  and  two  symmetric 
levels  V - V , V + V?  about  V_  as  shown  in  Figure  15a.  In  addition,  a fourth 
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level  V.,  can  also  be  defined  which  must  be  higher  than  + VT • The 
n U L 

amplitudes  V^,  and  as  well  as  the  polarity,  are  independently 
selectable.  The  filtered  Analog  INput  (AIN)  must  cross  from  below 

and  then  VQ  + from  below,  in  that  order,  before  the  Digital  INput  (DIN) 
is  raised.  A subsequent  crossing  of  from  above  will  cause  DIN  to  go  low. 
This  test  sequence  must  be  passed  by  the  analog  input  AIN  on  every  cycle . 
otherwise  the  burst  is  rejected,  the  logic  is  reset,  and  the  processor  waits 
for  the  next  burst.  If  the  analog  input  AIN  exceeds  V at  any  time  during 
the  burst,  the  burst  is  rejected,  the  logic  is  reset,  and  the  processor  waits 
for  the  next  burst  (cf . Figure  15a)  . Thus  the  amplitudes  Vq,  and 
specify  the  range  of  particle  sizes  accepted  by  the  processor  for  the  mea- 
surement. The  frequency  of  the  signal  DIN  can  be  prescaled  (cf . Figure  15b) 
to  produce  the  signal  DIN'  which  has  one  negative  slope  for  every  m negative 
slope  of  DIN.  The  constant  m is  front  panel  selectable  and  can  be  set  to 
1,  2,  4 or  8.  The  Analog  Input  Processor  is  designed  with  high-speed  ECL 
circuitry  and  stripline  technology  with  a maximum  Doppler  frequency  of  170  MHz. 

2 . The  Digital  Input  Processor 

The  prescaled  digital  input  DIN ' is  subsequently  processed  by  the  Digital 
Input  Processor.  The  negative  slope  of  DIN'  fires  a resectable  one-shot  T1 
of  duration  t^,  as  illustrated  in  Figure  16.  The  negative  slope  of  T1  fires 
a resettable  one-shot  of  duration  t^.  The  first  negative  slope  of  DIN'  also 
raises  a signal  G,  which  gates  a 100  MHz  crystal  clock  into  the  Flight  Time 
(FT)  counter,  as  illustrated  in  Figure  17.  The  next  negative  slope  of  DIN', 
retires  the  one-shot  Tl,  resets  the  T2  one-shot  and  generates  a short  TRansfer 
Pulse  (TRP)  which  stores  the  FT  buffer  and  increments  the  cycle  counter  (cf. 
Figures  16  and  17).  If  the  T2  one-shot  times  out  before  it  is  reset,  an  error 
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flag  Is  raised,  the  logic  is  reset  and  the  processor  waits  for  the  next 
burst.  As  a consequence,  the  logic  checks  that  every  period  of  the  pre- 
scaled digital  input  falls  within  predetermined  limits.  If  t is  the  period 
of  the  Doppler  signal,  we  must  have,  as  a consequence, 

< E1T  < + T2  , (43) 


where  m = 1,  2,  4,  8.  The  time  intervals  and  T ^ are  front  panel  selectable 
and  cover  the  range  60  nsec  < t^,  < 3 msec.  The  maximum  Doppler  frequency 

that  can  be  handled  by  the  Analog  and  Digital  Processors  working  in  tandem 
is  given  by. 


<V-> 


max(m) 

min(T^) 


8 

60  x 10 


- 130  MHz 


(44) 


3.  The  Timer 

The  100  MHz  crystal  clock  is  divided  by  10  to  form  a 10  MHz  time 

standard  which  is  used  as  the  clock  for  the  timer,  as  illustrated  in 

N 

Figure  17.  The  timer,  in  turn,  has  a 4 10  prescaler  which  sets  the 
time  base.  This  increments  a free  running  32-bit  counter  which  can  be 

a 

reset  if  required  but  is  normally  allowed  to  overflow.  The  gate  signal 
G is  inverted  to  form  G and  is  used  to  latch  the  contents  of  the  timer 
counter  into  the  timer  buffer,  as  illustrated  in  Figure  18.  The  purpose 
of  this  information  is  to  record  the  real  time  of  each  scattering  event, 
to  permit  the  reconstruction  of  the  velocity  fluctuations  in  time  (cf. 
Dimotakis  and  Brown  (1975)).  While  the  present  data  rate  was  often  high 
enough  (~  50  kHz)  to  warrant  an  examination  of  the  velocity  fluctuations 
in  time,  it  was  decided  not  to  Include  such  data  in  the  present  report. 
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4.  Performance  and  Operating  Modes 
The  LDV  processor  operates  at  the  theoretical  performance  of  r 1/2 
clock  count  (10  nsec)  in  determining  the  flight  time  &t^  of  the  i*''1 
scattering  particle.  The  overall  measurement  accuracy  is  of  course  limited 
by  a variety  of  other  factors  such  as  optical  precision,  signal- to-noise 
of  the  signal  at  the  processor  input,  etc.  In  all  the  high-speed  cases 
(M^  0 .8),  the  finite  frequency  of  the  processor  clock  was  the  limiting 
factor  in  the  velocity  measurement  accuracy  of  a single  particle.  The 
highest  frequency  that  occurred  was  27  MHz  in  the  Mes>*2.2  flow.  This  fre- 
quency was  measured  with  an  overall  accuracy  for  a single  particle  of  — 2%. 

A histogram  from  a 1,024  particle  record  taken  in  the  free  stream  (Me=“2.2) 
of  the  flight  time  for  8 Doppler  cycles,  is  depicted  in  Figure  19.  It  can 
be  seen  that  — 2/3  of  all  the  counts  fall  in  one  10  nsec  bin.  At  the  lower 
velocities  the  limitations  arose  as  a result  of  the  combined  effect  of  all 
the  other  sources  of  error.  For  the  0,1  measurements  an  accuracy  of 
0.5%  for  a single  particle  was  achieved.  It  should  be  emphasized  that  while 
the  accuracy  for  a single  particle  velocity  measurement  does  not  limit  the 
accuracy  with  which  mean  velocities  can  be  determined,  it  contaminates  the 
measurement  of  the  rms  fluctuation  levels.  *> 

The  processor  is  capable  of  operating  in  a variety  of  modes.  TVo  of 
the. modes  that  were  used  in  this  work  will  be  discussed. 

In  the  first  mode,  a fixed  number  of  cycles  is  stored  in  the  Mode  Logic 
circuitry  and  compared  with  the  contents  of  the  Cycle  Counter  after  each  TRP 
pulse  (cf . Figures  16  and  17) . When  the  two  numbers  match,  the  Count  End 
(CEND)  flag  is  raised  which  in  turn  lowers  the  GO  disabling  the  system 
(cf,  Figures  17  and  15a).  The  output  is  then  read  and  the  processor  is 
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reset  and  waits  for  the  next  burst  (dead  time  ~ 60  naec) . All  the  data 
for  ]^>  0.6  (JPL  20-inch  Wind  TUnnel)  were  recorded  in  this  mode.  Eight 
cycles  were  counted  for  the  beam  pair  at  0°  (u^  data)  with  the  prescaler 
m set  to  2,  Four  cycles  were  counted  for  the  ± 45°  (u  , u data)  beam 
pairs  with  the  prescaler  set  to  unity.  The  advantage  of  this  scheme  was 
that  the  digital  processor  saw  the  same  frequencies  and  numbers  of  cycles 
for  all  configurations.  This  occurs  for  the  present  measurements  because 
the  projection  of  the  velocity  vector  along  the  ± 45°  directions  is  approxi- 
mately given  by 


u 


2 


(45) 


and  hence  the  resulting  Doppler  frequencies  are  given  by 


1 

2 


(46) 


This  uniformity  in  the  processing,  allows  all  the  signals  to  be 
processed  with  identical  settings  except  for  a change  in  the  prescaler 
from  l.to  2.  The  particle  velocity  is  then  determined  by  reading  the 
output  of  the  FT  counter  which  corresponds  to  the  flight  time  of  the 
particle  for  the  fixed  number  of  fringes  specified. 

In  the  second  mode,  the  processor  requires  that  a fixed  number  of 
cycles  be  reached  or  exceeded  by  the  burst  and  records  the  flight  time 
and  number  of  cycles  of  the  entire  burst.  The  sequence  depicted  in 
Figure  15a  represents  such  a mode  where,  assuming  that  each  period  of 
DIN  passed  the  t^,  Tj  test,  the  end  of  the  burst  was  determined  by  the 
Analog  Input  Processor.  The  flight  time  At  that  would  have  resulted  from 
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such  a burst  is  indicated  on  the  figure  while  the  output  of  the  cycle 
counter  would  have  been  n * 3 (number  of  TRP  pulses  and  also  number  of 
fringe  plane  intervals).  Both  n and  At  would  be  recorded  in  this  mode 
and  the  velocity  would  be 


u 

x 


ns 

At 


(47) 


The  data  recorded  at  0.1  (GALCIT  Merrill  Wind  Tunnel)  were  acquired 
in  this  mode  requiring  at  least  four  fringes  for  the  ± 45  data  (u^,  u^) 
and  at  least  eight  for  the  0°  data  (u^) . Thus,  in  both  modes,  the  same 
minimum  number  of  fringes  had  to  be  crossed. 

There  are  several  important  differences  between  the  two  modes. 

First,  using  the  facility  of  the  processor  to  recognize  the  end  of  the 
burst  while  retaining  the  previous  valid  measurement  allows  the  total 
flight  time  to  be  used,  instead  of  an  arbitrary  fraction.  In  determining 
the  velocity  component  of  interest.  This  improves  the  measurement  accuracy 
in  an  obvious  way.  Second,  the  phase  errors  that  are  made  in  assigning 
equal  phase  to  the  crossings  of  the  reference  level  Vq  by  the  signal  are 
cancelled  when  the  entire  burst  is  used.  Third,  the  sampling  statistics 
are  different  between  the  two  modes.  In  the  first  mode  (fixed  number  of 
fringes),  all  we  know  is  that  the  particle  crossed  at  least  that  many.  In 
the  second  mode  (number  of  fringes  unrestricted  provided  it  exceeds  a certain 
minimum)  we  know  that  the  particle  crossed  exactly  that  many.  The  sampling 
bias  is  different  for  these  two  cases,  as  will  be  discussed  in  Section  VI-B. 

B.  Output  Formatting 

The  (i)  flight  time,  (ii)  total  number  of  validated  cycles  and  (iii) 
the  real  time  are  encoded  by  the  processor  into  three  32-bit  words  as  binary 
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coded  decimals  (BDC) . These  data  serve  as  input  to  the  subsequent  format- 
ting electronics  which  selected  the  four  least  significant  decades  of  the 
flight  time  to  form  a sixteen-bit  word  and,  depending  on  the  fringe  mode, 
either  the  sixteen  bits  of  the  real  time  for  the  fixed-fringe  mode,  or  the 
number  of  fringe  crossings,  in  the  free-fringe  mode.  Thes^e  data  were  then 
combined  to  form  a single  32-bit  word. 

The  asynchronous  data  from  the  processor  were  recorded  in  dual  1024  x 
32-bit  buffers.  One  buffer  was  filled  asynchronously  with  the  processor 
output  while  the  second  was  clocked  synchronously,  as  a sequence  of  4 bits/ 
word  on  a Kennedy  Model  9100  digital  tape  deck.  Each  1024  x 32-bit  buffer 
constituted  a 4K-bit  record  on  tape.  Several  records  were  recorded  at  each 
station  depending  on  the  turbulence  level  and  average  data  rate. 

VI.  Data  Processing 

In  computing  the  Reynolds  stress  -pu'v'  from  the  measured  velocity  data, 
it  is  necessary  to  estimate  both  the  first  and  second  moments  of  the  data, 
and  to  extract  relatively  small  differences  between  two  large  numbers.  This 
necessitates  a careful  accounting  of  all  random  and  systematic  errors  that 
result  from  spurious  measurements.  In  the  present  experiments,  such  errors 
are  caused  by  a finite  signal-to-noise  ratio,  and  appear  as  isolated  counts 
when  the  measurements  are  treated  in  histogram  form.  Errors  are  also  caused 
by  LDV  sampling  bias  and  by  both  the  finite  clock  period  and  by  noise  in  the 
processing  electronics  and  in  the  optics.  Errors  which  arise  as  a result  of 
the  finite  measurement  volume  have  been  found  to  only  be  important  in  the 
immediate  vicinity  of  the  wall,  and  have  not  been  corrected  in  the  present 
analysis.  The  correction  of  such  errors  would  result  in  a correction  of  at 
most  the  first  two  data  points  adjacent  to  the  wall,  for  which  we  have  the 
least  confidence,  and  would  not  affect  the  conclusions  presented. 
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A.  Histogram  Pruning 

The  data  from  the  present  experiments  were  recorded  as  multiple  records 
of  1024  individual  realisations  of  the  three  instantaneous  velocity  com- 
ponents u,  u + v,  u - v.  These  data  were  processed  by  forming  a separate 
histogram  for  each  record,  and  subsequently  performing  the  statistical 
computations  using  the  resulting  histograms.  For  the  high-speed  flow  experi- 
ments, the  data  were  recorded  in  the  fixed- fringe  mode  and  the  histograms 
were  formed  as  a function,  of  the  integer  flight  time  At^.  The  low-speed 
measurements  were  recorded  in  the  free-fringe  mode  for  which  the  flight 
time  was  large  and  covered  a much  wider  range.  For  these  data,  the  histogram 
was  formed  as  a function  of  the  velocity  of  the  particles,  computed  as  an 
integer  percentage  of  the  maximum  velocity  that  the  processor  would  admit, 
i.e . , 


I.  - 100  t1  n.  , (48) 

“i 

where  is  the  minimum  Doppler  period  and  n^  is  the  number  of  fringes 
crossed  by  the  iCh  particle  in  time  At^. 

The  pruning  was  performed  in  two  passes.  First,  data  in  any  bin  that 
contained  only  one  count,  ar.d  did  not  have  neighboring  bins  with  more  than 
one  count,  were  discarded.  In  the  second  pass,  any  data  isolated  from  the 
main  body  of  the  histogram  by  at  least  one  zero  were  also  discarded.  If, 
as  a result  of  these  two  operations,  more  than  24  measurements  out  of  the 
total  of  1024  were  rejected,  the  entire  record  was  rejected.  This  scheme 
was  preferred  to  the  more  usual  one  of  rejecting  data  outside  a certain 
number  of  standard  deviations  from  the  mean.  By  way  of  example,  the  few 
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counts  in  bin  40  of  the  histogram  in  Figure  19  would  be  discarded  by  this 
scheme. 

The  histograms  which  result  from  the  data  in  the  present  investigations 
are  illustrated  in  Figures  20  and  21.  In  these  figures,  the  histograms  at 
several  locations  through  the  boundary  layer  are  superimposed  on  the  mean- 
velocity  profiles  for  that  case.  For  both  flows  represented  here,  the 
histogram  for  the  free-stream  flow  encompasses  only  four  bins.  Hear  the 
wall,  the  histogram  is  spread  over  more  than  forty  bins. 

B.  Sampling  Bias 

The  complexity  of  processing  single  particle  laser-Doppler  velocity 
data  is  compounded  by  the  fact  that  the  fluid  velocity  is  sampled  in  a 
biased  manner.  This  problem,  first  pointed  out  by  McLaughlin  and  Tiederman 
(1973),  is  a consequence  of  the  fact  that  the  particle  flux  through  the 
measurement  volume,  and  the  resulting  measurement  probability  per  unit  time, 
is  higher  when  the  local  fluid  velocity  is  high  than  when  it  is  low.  The 
fact  that  the  measurement  probability  is  a function  of  the  measured  quantity 
gives  rise  to  the  sampling  bias  p(u),  where  u is  the  three-dimensional  velocity 
vector.  Consequently,  to  compute  the  expectation  value  of  a particular 
function  of  the  measurements  we  must  remove  this  bias,  i.e,. 


(49) 


The  dependence  of  the  sampling  bias  function  P on  the  local  velocity  vector 
u,  on  the  shape  of  the  measurement  volume,  and  on  the  minimum  number  of 
fringe  crossings  required  by  the  processing  electronics,  has  been  derived 
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by  Dimotakis  (1976).  For  data  recorded  in  the  fixed-fringe  mode,  if  the 
angle  between  Che  intersecting  beams  is  sirall,  the  expression  for  0(u) 
simplifies  to 

1./2 

P(uje)  “ u^(l  + u^/u^2)  jjl  - e2(l  + u„2/u^  )]  i (50) 

where  is  the  velocity  component  perpendicular  to  the  fringe  planes, 
is  the  velocity  component  parallel  to  the  fringe  planes  and  perpen- 
dicular to  the  beam  bisector,  and  e is  the  ratio  n . /n„  of  the  minimum 
number  of  fringe  crossings  required  by  the  processor  to  the  total  number 
in  the  focal  volume.  Note  that  recording  the  flight  time  for  a fixed 
number  of  fringes  in  fact  requires  that  at  least  that  many  were  crossed. 
Consequently,  the  fixed  number  of  fringes  is  equal  to  n ^ in  the  context 
in  which  e was  defined  for  Equation  (50). 

Evaluating  0(u;e)  for  the  data  in  the  boundary  layer,  which  were 

acquired  using  the  fixed-fringe  mode  (0.6  ~ 2,2),  we  note  that 

,2  2 2 

u * /u  = v /u  « 1 so  that,  in  this  case, 
a 1 

P(u;e)  -uA'(l  - e2)  , (51) 

or, 

0(u;e)  « ux  ■ (52) 

In  going  from  Equation  (51)  to  Equation  (52),  an  additional  assump- 
tion has  been  made.  Even  though  n is  fixed  by  the  signal  processing 
electronics,  the  total  number  of  fringes  crossed,  n^,,  is  unfortunately  a 
complicated  function  of  the  particle  size.  This  is  a consequence  of  the 
fact  that  the  amplitude  levels  VQ  and  are  absolute.  A large  particle. 
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with  a correspondingly  large  scattering  cross-section,  will  yield  a burst 
which  is  expected  to  result  in  more  validated  cycles.  Thus,  the  ratio  e 
is  a function  of  the  particle  size.  In  the  context  of  Equation  (49),  the 
transition  from  Equation  (51)  to  (52)  is  valid  only  if  the  local  velocity 
vector  is  uncorrelated  with  particle  size.  This  requirement  will  be  satisfied 
if  the  particles  track  the  flow. 

With  this  proviso  we  then  have  from  Equations  (49)  and  (52),  for  the 
u ^ data  in  the  fixed- fringe  mode, 


<f(u)>  =■ 


(53) 


In  particular,  the  mean  velocity  is  given  by 


<ux>  = 


(54) 


where  N is  the  total  number  of  data  points,  as  originally  proposed  by 
McLaughlin  and  Tiederman  (1973).  Similarly  the  mean  square  is  given 


i 


(55) 


The  mean-square  fluctuations  can  then  be  computed  by  means  of 

<u'x2>  = <ux2>  - <ux>2  . 


(56) 
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1 


The  latter  form  is  preferable  in  avoiding  finite  precision  problems  in 
the  computer,  at  the  expense,  however,  of  substantially  more  computing 
time. 

In  evaluating  Equations  (54)  and  (55),  note  that  the  i^  measured 
velocity  component  is  given  by  Equation  (47) , or 


u , . = 


(58) 


where  n is  equal  to  the  number  of  fringes,  s is  the  fringe  spacing,  is 
the  clock  period  (10  nsec),  k^  is  the  (integer)  output  of  the  flight  time 
counter  and  x is  a random  variable  representing  the  collective  uncertainty 
in  the  determination  of  the  flight  time  at  and  resulting  from  (i)  the 
finite  clock  period  , and  (ii)  the  uncertainty  with  which  the  equal 
phase  points  are  determined  by  the'  comparators  from  the  signal  burst  as 
a result  of  finite  s ignal-to-noise  ratio.  The  1/2  is  added  to  k^  to  remove 
the  systematic  truncation  error  of  the  digital  counter. 

Using  Equations  (54)  and  (55)  and  correcting  for  the  finite  variance 
of  the  random  variable  x,  we  have 


and 


/ v _ / ns  \ K 

'“infixed  " (^Tcj  J_L 

■(?)' 


(59) 


(60) 
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where 


* 


and 


J x2p(x)dx 
J p(x)dx 


(61) 


(62) 


is  the  variance  of  x. 

Unfortunately,  the  situation  is  less  well  defined  for  the  ± 45°  data 
recorded  in  the  fixed-fringe  mode.  In  that  case,  the  two  relevant  veloci- 
ties, ux  and  u,,  are  of  the  same  order  of  magnitude.  Nevertheless,  since 
only  one  velocity  component  was  measured  at  a time,  we  use  Equation  (52) 
for  the  bias  and  reduce  the  ± 45°  data  using  Equations  (59)  and  (60) . While 
this  is  clearly  not  correct,  it  is  better  than  ignoring  the  bias  altogether. 

Note  that  the  expression  for  the  sampling  bias  0 given  by  Equation  (52) 
is  not  appropriate  for  data  acquired  in  the  free-fringe  mode  because  this 
equation  was  derived  using  the  assumption  that  the  probability  of  measure- 
ment is  determined  by  the  fact  that  the  particle  has  crossed  at  least  the 
fixed  number  of  fringes.  In  the  free-fringe  mode,  the  measurement  of  the 
number  of  fringes  crossed  and  of  the  flight  time  in  units  of  t corresponds 
to  a particle  crossing  of  exactly  fringes.  The  probability  for  this  event 
can  be  determined  from  Equation  (50-) , Since  |3  (u;  n^/n^,)  is  proportional  to 
the  probability  of  crossing  at  least  n^  fringes,  i.e., 

&fixed  «P  {n<“il“l  » (&3) 
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we  have  that 


[n  = ^1  u}  = p [n  < n±  | u]  - p [n  < {n1  + 1)  j u]  , 


(64) 


where  p { } denotes  the  probability  per  unit  time  of  observing  the  event 
in  the  braces.  Consequently, 

i 

^free^3  W “ Pfixed(^;ni/V  ' 3fixed[^(ni  + * (63) 


or 


3/2 

Pfree*-;,W  * UJ.  + U/<Z>  (ni  + 1/1 > 


(66) 


where,  in  the  context  of  Equation  (49),  the  assumption  has  been  made  that 
the  particle  size  is  uncorrelated  with  the  local  velocity.  In  the  context 
of  Equation  (49),  note  that 


(67) 


whereas 


<ui  + u/) 


> 


(68) 


where  At^  is  the  time  of  flight  for  the  measured  n^  fringes.  Substituting 
these  equations  in  Equation  (66) , we  have 


-1 


3free  (^ni/nT>  06 


(k^  + 1/2  + x)n^ 
ni  + 1/2 


(69) 


where,  as  before,  ^ is  the  integer  output  of  the  flight  time  counter  and 
x is  a random  variable  representing  the  uncertainty  in  the  measurement  of 
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the  flight  time.  Using  this  bias  function,  we  thus  have 


l v _ /na\  Vj, 

1>f"*  v«J  K-i,: 


(70) 


and 


(71) 


where 


(^i  + j)^^  + 


(72) 


and  cT^  is  the  variance  of  x,  Note  that  Equations  (70)  and  (71)  are  equally  ^ 
valid  for  the  0°  data  as  for  the  ± 45°  data.  The  data  recorded  at  0‘1- 
in  the  free-fringe  mode  were  analyzed  using  Equations  (70)  and  (71).  In 
order  to  compare  them  to  the  previous  data,  Equations  (54)  and  (55)  were 
also  used,  which  in  this  case  reduce  to 


<v 


(73) 


and 


(74) 


where 


L. 

£,m 


m 

n . 

1 


<ki + V 


(75) 
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Note  also  that  Equation  (69)  for  the  sampling  bias  in  the  free-fringe 

mode  is  in  disagreement  with  the  results  of  Hoesel  and  Rodi  (1977)  and 

Buchhave,  George  and  Lumley  (1979),  who  claim  that  5 oc  &t^  « (k^  + ^) . 

If  that  were  true.  Equations  (59)  and  (60)  would  be  valid  for  the  free- 

fringe  mode  as  well,  with  the  three  summations  as  defined  in  Equation  (61). 

2 

The  variance  c of  x was  estimated  from  the  mean-square  fluctuations 

2 

in  the  free  stream  where  (u  ' ) ^ 0. 

The  results  of  this  analysis  have  been  applied  to  the  data  obtained 
in  these  experiments.  The  bias  correction  is  most  important  in  the  region 
near  the  wall,  where  the  fluctuation  levels  are  highest.  These  calculations 
indicate  at  most  a 2%  correction  for  the  low-speed  flow  data  in  this  region, 
and  a 3.3%  correction  for  the  high-speed  flow  data. 

VII . Results  and  Discussion 
A.  Data  Reduction 

Following  the  computation  of  the  proper  statistical  averages,  the  data 
were  reduced  to  a form  appropriate  for  boundary  layer  flows  by  utilizing  the 
scheme  outlined  in  Part  I of  this  report  by  Collins,  Coles  and  Hicks  (1978). 

For  this  purpose,  the  Van  Driest  (1955)  scaled  mean  streanwise  velocity  data 
were  fitted  to  the  law  of  the  wall  and  the  law  of  the  wake  using  a three- 
parameter  nonlinear  least-squares  fit  with  a single  constraint  equation  based 
on  the  local  friction  law.  The  profile  parameters  u^,  II  and  6,  determined 
by  this  analysis,  were  then  used  to  infer  both  the  distribution  of  shear  stress, 

I 

and  of  the  normal  velocity,  from  the  data  by  an  integration  of  the  equations 
for  the  conservation  of  mass  and  momentum  from  the  wall  to  the  free  stream, 
using  the  formulation  by  Spalding  (1961)  to  describe  the  flew  in  the  viscous 
sublayer.  The  details  of  these  computations  are  given  in  Part  1. 
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One  significant  difference  between  the  present  LDV  data  and  the  data 
obtained  using  standard  Pitot  instrumentation  (cf.  Part  I)  is  the  ability 
of  the  Pitot  to  obtain  valid  data  in  the  region  adjacent  to  the  wall. 

Because  of  noise  considerations , the  present  data  are  limited  to  a distance 
of  approximately  one  focal-volume  diameter  from  the  wall.  This  fact  requires 
that  the  boundary- layer  scale  parameters  6*  and  0 be  computed  directly  from 
the  law  of  the  wall  and  the  law  of  the  wake  representations  rather  than  from 
a direct  integration  of  the  data,  as  in  the  results  discussed  in  Part  I. 

The  scale  parameters  are  defined  as 


and 


(76) 


(77) 


In  terms  of  the  Van  Driest  (1955)  scaled  velocity  U, 
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and  U is  defined  by  the  Van  Driest  scaling 
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Thus, 
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The  integrals  involved  in  these  expressions  were  evaluated  in  terms 
of  the  integrals  P and  Q defined  in  Part  I.  These  integrals  are  functions 
only  of  the  mean-velocity  profile  parameters  u^,  E and  6,  and  are  evaluated 
in  terms  of  the  Spalding  formulation  for  the  boundary  layer. 

As  a result  of  the  differences  in  the  computational  procedure  for  the 
scale  parameters  6*  and  S,  these  quantities  differ  slightly  from  those 
computed  for  the  identical  flows  given  in  Part  I of  this  work.  Corrections 
for  these  discrepancies  have  been  included  where  comparisons  are  made  to 
the  previous  work. 

B.  Streamwise  Mean-Velocity  Profiles 

A comparison  between  the  mean-velocity  profiles  obtained  in  the  present 
experiments  and  those  obtained  from  the  Pitot  data  of  Part  I is  given  in 
Figures  22  and  23.  In  order  to  make  a consistent  comparison,  the  Pitot  data 
from  Part  I have  been  replotted  using  the  integral  scale  0 obtained  from  the 
LDV  experiments.  At  high  Mach  numbers,  the  results  exhibit  a small  discrep- 
ancy between  the  LDV  and  the  Pitot  data,  with  the  LDV  data  exhibiting  lower 
mean  velocities  near  the  wall.  The  cause  of  this  discrepancy  is  not  clear 
at  this  writing,  nor  is  it  clear  which  measurement  is  the  more  reliable. 

The  low-speed  flow  mean-velocity  profiles  have  been  computed  using  both  the 
exact  bias  correction,  according  to  Equation  (70)  and  using  the  harmonic 
mean  according  to  Equation  (59).  The  differences  are  small  as  are  the 
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differences  between  these  results  and  the  results  of  the  mean  value  of 
these  data  without  any  bias  correction.  The  comparisons  have  not  been 
plotted  here.  The  agreement  between  the  LDV  data  and  the  Pitot  data  is 
excellent,  for  the  low-speed  results. 

The  data  obtained  from  the  present  experiments  have  been  plotted 
in  coordinates  appropriate  to  the  law  of  the  wall  in  Figures  24  and  25. 

The  solid  line  associated  with  each  data  set  has  been  computed  using  the 
formulation  by  Spalding  (1961),  and  illustrates  the  fit  obtained  in  the 
computation  of  u^,  !1  and  6.  The  high-speed  flow  results  shown  in  these 
figures  illustrate  the  problem  of  obtaining  measurements  near  the  wall. 
These  data  exhibit  a departure  from  the  computed  profile  because  of  noise 
which  occurs  when  the  focal  volume  intercepts  the  wall.  This  fact  required 
that  the  data  included  in  the  least-squares  fit  be  limited  to  y 5 200  for 
these  data,  as  was  the  case  for  the  Pitot  data  reported  in  Part  I. 

The  profile  parameters  and  5 computed  from  the  least-squares 

procedure  agree  closely  with  those  parameters  computed  from  the  Pitot  data 
for  the  same  flows.  The  computations  of  the  wake  parameter  II,  however,  are 
considerably  higher  in  every  case.  The  self-consistency  of  the  LDV  results 
indicates  that  this  discrepancy  is  probably  not  the  consequence  of  data 
reduction  errors . 

C.  Streanwise  Velocity  Fluctuations 

The  fluctuations  of  the  streairwise  velocity  component,  normalized  by 
the  free-stream  velocity,  are  shown  in  Figures  26  and  27  as  a function  of 
the  normal  coordinate.  These  results  are  compared  to  the  distribution  of 
the  u/  fluctuations  obtained  at  low  speeds  by  Klebanoff  (1954)  in  flow  at 
constant  pressure,  represented  by  the  solid  squares  in  these  figures.  For 
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these  comparisons,  the  normal  coordinate  for  the  data  by  Klebanoff  have 
been  scaled  with  the  boundary- layer  scale  obtained  in  the  present  experi- 
ments . 

The  error  bars  in  each  of  the  subsequent  figures  have  been  computed 
as  the  standard  deviation  of  the  ensemble  of  records  used  in  computing 
the  mean  values.  In  each  case,  a few  data  points  for  which  the  error 
bars  exceed  the  mean  value  have  been  excluded.  From  these  results,  it 
is  clear  that,  near  the  tree-stream,  the  velocity  fluctuations  are  properly 
computed  from  the  LDV  data.  As  the  wall  is  approached,  the  deviations 
between  the  present  results  and  those  of  Kelbanoff  are  within  the  confidence 
limits  in  the  present  data. 

It  can  be  concluded  from  these  data  that,  within  the  error  in  the 
present  experiments,  no  measurable  effects  of  compressibility  have  been 
found  on  the  mean-square  velocity  fluctuations.  The  data  for  Mach  numbers 
up  to  2.2  are  adequately  represented  by  the  data  of  Klebanoff.  The  differences 
exhibited  in  these  figures  between  the  high-speed  flow  data  and  the  low-speed 
flow  data  arise  from  the  differences  in  the  boundary- layer  Reynolds  number, 

Re^.  These  results  confirm  Morkovin's  (1961)  hypothesis  that  there  is  no 
essential  difference  in  the  dynamic  behavior  of  the  boundary  layer  at  constant 
pressure  for  Mach  numbers  up  to  4 or  5. 

D.  Reynolds  Stress  Profiles 

The  distribution  of  Reynolds  stress,  computed  as  - pu  'v ' , is  shown  in 
Figures  28  and  29  for  the  data  of  these  experiments.  In  order  to  compare 
the  Reynolds  stress  deduced  from  the  experimental  data  with  the  expected 
distribution  through  the  boundary  layer,  a second  computation  based  on  the 
integral  formulation  outlined  in  Part  I is  displayed  as  the  solid  curve  in 
each  figure.  The  results  from  the  integral  formulation  coincide  with  the 
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results  obtained  using  the  Pitot  data,  with  minor  differences  arising 

V 

from  the  difference  in  the  computed  value  of  the  wake  parameter  II. 

For  comparison  with  the  expected  behavior  at  low  speeds,  the  data 
of  Klebanoff  (1954)  has  been  Included  in  these  figures  as  the  dark  squares 
on  each  figure.  The  agreement  between  the  integral  formulation  and  the 
data  of  Klebanoff  indicates  that  the  integral  formulation,  based  on  the 
mean-velocity  data,  yields  the  expected  distribution  for  the  total  shear 
stress . 

The  dominant  feature  of  the  results  for  the  distribution  of  the 
Reynolds  stress,  as  measured  by  the  laser-Doppler  technique,  is  the  de- 
parture of  the  measured  Reynolds  stress  from  the  results  of  the  integral 
computation  in  the  region  y/8  ~ 3.  This  departure  represents  a major  error 
in  these  measurements  and  may  indicate  a fundamental  limitation  on  the 
applicability  of  the  laser-Doppler  technique  in  flows  of  this  type. 

Several  conclusions  can  be  reached  from  these  data.  The  anomaly 
represented  by  the  departure  of  the  Reynolds  stress  from  the  expected  value 
near  the  wall  has  been  reported  previously  by  Johnson  and  Rose  (1973),  by 
Yanta  and  Lee  (1974),  and  by  Abbiss  (1976).  An  explanation  of  this  phe- 
nomenon by  Sandborn  (1974)  assumed  that  the  phenomenon  is  related  to  the 
contributions  of  density  fluctuations  to  the  Reynolds  stress.  That  this 
assumption  is  incorrect  is  clear  from  the  present  data.  The  observed  phe- 
nomenon is  not  Ifech  number  dependent  and  hence  is  not  an  effect  of  com- 
pressibility. These  data  also  indicate  that  the  discrepancy  in  the  Reynolds 
stress  is  not  associated  with  particle  lag  at  the  position  of  measurement, 
because  for  this  to  be  the  case  the  results  of  Equation  (13)  require  that 
the  lag  occur  at  increasing  distances  from  the  wall  with  increasing  free- 
stream  velocity,  In  addition,  the  agreement  obtained  in  the  distribution 
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of  the  mean-square  velocity  fluctuations  indicates  that  the  particles 
follow  the  flow  locally  throughout  the  boundary  layer. 

Similar  distributions  for  the  Reynolds  stress,  measured  using  hot- 
wire anemometer  techniques  by  Laderman  (1978),  have  recently  been  shown 
to  agree  with  the  present  results  given  by  the  integral  distributions 
shown  by  the  solid  curves  in  Figures  28  and  29. 

These  results  again  confirm  Morkovin1 s hypothesis  (1961)  and  refute 
the  hypothesis  by  Sandborn  (1974)  that  density  fluctuations  may  be  re- 
sponsible for  major  changes  in  the  distribution  of  Reynolds  stress  at 
high  speeds.  Similar  conclusions  were  expressed  in  Parc  I,  based  on 
the  results  of  the  Pitot  measurements  of  the  mean-velocity  profiles. 

E,  The  Problem  of  the  Normal  Velocity 

The  results  for  the  distribution  of  normal  velocity  are  illustrated 
in  Figure  30  for  the  high-speed  flow  at  low  Reynolds  number.  In  this 
figure,  the  solid  line  is  the  result  for  the  correct  distribution  of 
v/u,  computed  from  the  integral  representation  of  the  data  outlined  by 
Collins,  Coles  and  Hicks  (1978),  using  the  boundary- layer  parameters  uT , 

II  and  6 derived  from  the  least-squares  analysis. 

The  errors  in  estimating  the  normal  velocity  from  the  present  data 
are  substantial  in  every  case.  The  principal  feature  of  the  data  is  that 
the  normal  velocity  is  large  and  negative  at  the  wall,  increasing  sharply 
in  magnitude  as  the  wall  is  approached.  The  disparity  between  the  large 
negative  value  and  the  error  bars  indicates  that  the  particles  are  being 
convected  toward  the  wall  in  this  region.  In  the  free  stream,  the  velocity 
is  positive  in  most  cases,  but  with  a magnitude  that  is  larger  than  the 
correct  free-stream  value  as  represented  by  the  integral  formulation. 
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These  data  illustrate  the  difficulty  in  accurately  representing 
the  normal  velocity.  This  difficulty  arises  as  a consequence  of  the 
fact  that  any  small  misalignment  of  the  optics  will  result  in  substantial 
errors  in  the  computed  value  of  the  normal  velocity,  which  has  a maximum 
value 

HC. 

(v/u)  ~ 2 « 

where  H = 6*/6  is  the  boundary  layer  form  parameter.  In  addition,  because 
of  the  non-uniform  particle  distributions,  the  normal  velocity  components 
of  both  the  flow  near  the  wall  and  the  flow  in  the  vicinity  of  the  boundary- 
layer  edge  are  improperly  represented. 

The  fluctuations  in  the  normal  velocity  are  represented  for  these  same 
data  in  Figure  31.  These  data  have  approximately  the  value  of  the  data  of 
Klebanoff  (1954)  for  incompressible  flow,  represented  by  the  solid  squares 
in  the  figure.  However,  in  light  of  the  difficulties  in  measuring  the 
normal  velocity,  this  agreement  is  fortuitous.  The  magnitude  of  the  error 

c* 

bars  for  these  data  are  approximately  one-half  of  the  magnitude  of  the  data, 
and  indicate  a low  degree  of  confidence  in  the  results. 

VIII.  Conclus ions 

The  present  experiments  lead  to  a number  of  conclusions  regarding  the 
applicability  of  the  laser-Doppler  technique  to  the  study  of  boundary- layer 
flows.  The  results  presented  for  the  behavior  of  the  streanwise  velocity 
components  indicate  the  importance  of  including  in  the  analysis  of  the  raw 
data  an  adequate  technique  for  pruning  the  initial  histograms  in  order  to 
remove  errors  in  the  data  introduced  by  noise.  These  results  further  indicate 
the  importance  of  including  proper  statistical  averaging  in  the  analysis  in 
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order  to  account  for  the  bias  introduced  by  sampling.  In  addition,  the 
effects  of  a finite  clock  period  and  noise  have  been  shown  to  be  Important 
in  the  proper  determination  of  the  streamwise  velocity  fluctuations,  uy, 
near  the  boundary-layer  edge.  When  these  effects  are  properly  included 
in  the  data  analysis,  a direct  comparison  between  the  present  laser-Doppler 
measurements  and  measurements  in  the  same  flows,  presented  in  Part  1 of 
this  work,  indicates  good  agreement  for  the  streamwise  velocity  components 
u and  u ,.  The  present  results,  when  properly  scaled  to  account  for  com- 
pressibility using  the  ideas  of  Van  Driest  (1955),  are  also  shown  to  agree 
with  the  low-speed  data  presented  by  KLebanoff  (1954),  This  observation 
further  corroborates  the  conclusions  expressed  in  Part  I of  this  work 
regarding  the  proper  role  of  compressibility  in  determining  the  behavior 
of  the  mean  velocity,  and  confirms  Morkovin's  hypothesis  (1961)  regarding 
the  role  of  density  fluctuations  in  the  boundary  layer  at  constant  pressure. 

The  present  measurements  of  the  Reynolds  stress,  - pu  'v  ,,  and  for  the 
normal  velocity,  v,  are  in  substantial  disagreement  with  the  expected  results 
derived  from  the  mean  flow.  These  observations  seem  to  indicate  a funda- 
mental difficulty  with  the  application  of  the  laser-Doppler  technique  as 
implemented  in  these  experiments  to  the  measurement  of  the  detailed  structure 
of  flows  in  the  vicinity  of  a wall.  The  measured  Reynolds  stress  profiles, 
-pu'v',  deviate  from  the  expected  behavior  for  y/0  3,  indicating  an  error 

in  the  measurement  for  the  u “v  ' correlation  in  the  region  adjacent  to  the 
wall.  This  error  arises  because  of  the  strong  correlation  between  the  normal 
velocity  component  and  the  particle  number  density  that  occurs  as  a result 
of  the  depletion  of  particles  near  the  wall.  These  results  indicate  that 
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the  flaw  Inside  the  viscous  sublayer  is  such  that  a particle  that  enters 

the  sublayer  has  a very  low  probability  of  leaving.  The  scaling  law  for 

particle  behavior  in  a turbulent  boundary  layer  (cf.  Equation, (13))  is  not 

applicable  in  the  viscous  sublayer.  A particle  which  enters  the  viscous 

sublayer  is1  subjected  to  the  very  high  Lagrangian  frequencies  which  occur 

in  the  motion  of  the  longitudinal  vorticies,  whose  transverse  extent  is 

approximately  , and  lagp  the  flow  in  this  region.  As  a consequence, 

UT 

a fluid  element  coming  from  the  wall  is  less  likely  to  carry  particles  than 
a fluid  element  moving  toward  the  wall.  This  results  in  an  under- 
representation of  the  positive  component  of  the  normal  velocity,  v,  in 
the  region  near  the  wall,  resulting  in  a negative  value  for  the  measured 
component  normal  to  the  wall. 

This  under-representation  of  the  motion  near  the  wall  also  results 
in  a lower  measured  value  for  the  Reynolds  stress,  -pu/v,J  than  would  be 
expected  by  an  analysis  based  on  the  mean  flow  as  described  in  Part  I,  since 
the  upswelling  of  fluid  from  the  wall  is  associated  with  a lower  streanwise 
velocity.  This  bursting  motion  has  been  observed  by  Blackwelder  and  Kaplan 
(1976)  to  be  intermittent  and  quite  violent,  and  is  held  to  be  responsible 
for  a large  fraction  of  the  total  stress  near  the  wall  ~ 100^.  The 

particle  transport  to  the  wall  via  the  viscous  sublayer  is  so  effective, 
even  at  the  lower  velocities  of  the  Merrill  wind  tunnel  experiments,  that 
the  wall  appears  as  an  infinite  sink  for  particles. 

We  conclude  from  these  observations  that  the  measured  disparity  between 
the  Reynolds  stress,  -pu'v',  and  the  total  stress,  as  defined  by  the  inte- 
gration of  the  mean  flow,  arise  not  as  a consequence  of  compressibility, 
as  suggested  by  Sandborn  (1974),  but  as  a consequence  of  particle  depletion 
which  is  unique  to  boundary  layers  in  air.  This  is  a fundamental  limitation 
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of  the  techniques  for  particle  seeding  which  must  be  solved  if  this 
technique  is  to  be  useful  for  detailed  measurements  near  the  wall. 

For  free-shear  flows,  because  of  the  absence  of  the  wall,  this  problem 
does  not  arise,  and  the  laser-Doppler  technique  will  give  an  accurate 
representation  of  the  mean  and  fluctuating  velocity  components  and  their 
correlations . 
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Figure  1.  The  Doppler  frequency  as  a function 
of  flow  velocity. 
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Figure  2.  The  number  of  fringes  in  the  focal  volume  as  a 
function  of  the  focal  volume  diameter. 
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Figure  4.  The  structure  of  the  turbulent  boundary  layer  as  visualized 
by  the  use  of  aluminum  flakes  in  water;  photograph  courtesy 
of  Brian  Cantwell,  Caltech  (1978). 
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Figure  5.  The  ratio  of  the  boundary- layer  edge  velocity  to 
the  boundary- layer  thickness  in  wall  coordinates 
as  a function  of  Re^, 
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Figure  6.  The  particle  size  distribution. 
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Figure  7.  The  distribution  of  particle  number  density, 
Refl  = 23,000. 
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Figure  8.  The  distribution  of  particle  number  density. 
Re  = 40,000. 


62 


1.  *.  3.  4.  B.  8*  7.  8*  8.  10*  1J . 12-  18.  H.  16 

Y/THETR 

Figure  9.  The  distribution  of  particle  number  density, 
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Figure  10 


High-speed  laser-Doppler  velocimeter 
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Figure  13.  Laser-Doppler  anemometer  signals  at  the  outer  edge  of  a 
turbulent  boundary  layer,  = 2.2,  Re^  = 40,000. 
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Figure  14.  LDV  processor  block  diagram. 
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Figure  16.  The  logic  associated  with  the  digital  Input  processor- 
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Figure  19,  A typical  histogram  from  the  free-stream 
flow  at  M - 2, 2. 
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Figure  20.  Representative  histograms  from  the  high- 
speed boundary  Layer. 
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Figure  23,  Mean  velocity  profiles. 

Re.  - 40,000. 
cJ 


77 


0 soefroog 


AEDC-TR-79-49 


n=2.lG7 


Figure  25.  Mean  velocity  profiles  with  Van  Driest  scaling, 

Re.  = 40,000. 
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Figure  29.  The  distribution  of  Reynolds  stress, 
Re0  ■ 40,000. 
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Nomenclature 


Symbol 

Equation 

Meaning 

c 

{22) 

constant  in  walL  law  (5.0) 

cf 

(84) 

local  friction  coefficient 

d 

(7) 

particle  diameter 

P 

k. 

integer  output  of  the  flight  time 

1 

counter 

n 

(21) 

particle  number  density 

n. 

l 

number  of  observations  with  velocity  u^ 

Re0 

Reynolds  number  based  on  9 

s 

(2) 

fringe  spacing 

U,  V 

strearawise  and  normal  velocity  component 

u . 

1 

(29) 

velocity  component  in  i^1  direction 

UT 

(18) 

friction  velocity 

u 

(1) 

velocity  component  normal  to  the 

L 

■ fringe  plane 

u 

velocity  component  parallel  to  the 

H 

fringe  plane 

X,  y 

, streanwise  and  normal  coordinates 

a, 

y,  0 Figure  12 

angles  defined  by  the  beam  planes 

P(u) 

(49) 

bias  function 

b 

(26) 

boundary- layer  thickness 

e 

(77) 

boundary- layer  momentum  thickness 

3i 

Figure  12 

beam  separation  angles 

K 

(22) 

Karman  constant  (0.41) 

X 

laser  wave  length  (0.5145  um) 

viscosity 

V 

kinematic  viscosity 
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Nomenclature  (Cont.) 


Symbol 

Equation 

Meaning 

VD 

(1) 

'Doppler  frequency 

n 

(22) 

strength  of  the  wake  component 

T 

C 

(5) 

processor  clock  period 

T 

P 

(8) 

particle  response  time 

T 

u 

shear  stress  at  the  wall 

Subscripts 


< >e 

edge  or  external  value 

< >f 

pertaining  to  the  fluid 

C>P 

pertaining  to  the  particles 

< >w 

wall  value 

( )' 

rms  fluctuations 

Superscripts 

( ) + 

value  made  dimensionless  with 

87 


